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ABSTRACT 

CX Cep (WR 151) is the WR+0 binary (WN5+05V) with the second short- 
est period known in our Galaxy. To examine the circumstellar matter distribution 
and to better constraint the orbital parameters and mass-loss rate of the WR star, 
we obtained broadband and multi-band (i.e. UBVRI) linear polarization obser- 
vations of the system. Our analysis of the phase-locked polarimetric modulation 
confirms the high orbital inclination of the system (i.e. i = 65°). Using the 
orbital solution of Lewis et al. (1993) we obtain masses of 33.9M0 and 23.9M0 
for the O and WR stars respectively, which agree with their spectral types. A 
simple polarimetric model accounting for finite stellar size effects allowed us to 
derive a mass-loss rate for the WR star of 0.3 — 0.5 x 10~^MQ/yr. This result was 
remarkably independent of the model's input parameters and favors an earlier 
spectral type for the WR component (i.e. WN4). Finally, using our multi-band 
observations, we fitted and subtracted from our data the interstellar polarization. 
The resulting constant intrinsic polarization of 3 — 4% is misaligned in relation 
to the orbital plane (i.e. 6cip = 26" vs. Q = 75°) and is the highest intrinsic 
polarization ever observed for a WR star. This misalignment points towards a ro- 
tational (or magnetic) origin for the asymmetry and contradicts the most recent 
evolutionary models for massive stars (Meynet & Maeder 2003) which predict 
spherically symmetric winds during the WR phase (i.e. CIP = 0%). 

Subject headings: Polarization - binaries: eclipsing - stars: Wolf-Rayet - tech- 
niques: polarimetric - stars: individual: CX Cep - stars: mass loss 
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1. Introduction 

Wolf-Rayet stars have always puzzled astronomers in their attempt to explain their 
extreme mass-loss and dense winds. Binarity was initially believed to be the only mechanism 
capable of expelling their outer layers, expose their CNO-enriched cores and produce their 
strong winds. However, the existence of single Wolf-Rayet stars with similar characteristics 
as their double-star counterparts was a big drawback for this scenario. 

Radiative (line-driven) wind acceleration (Cassinelli & Castor 1973) is now the widely 
accepted explanation for the Wolf-Rayet phenomenon in massive stars. The strong UV 
radiation produced by massive stars coupled with significant metallicity provides a suit- 
able environment for WR-type winds to develop. The inclusion of rotation in massive star 
evolutionary models is now able to reconcile theoretical predictions with previously unex- 
plained observational evidence (Maeder & Meynet 2000). For instance, the observed ratio of 
nitrogen-rich (WN-type) WR stars to carbon-rich (WC-type) WR stars is better reproduced 
by models that include rotation. 

One important prediction of stellar evolution with rotation is that massive stars with 
solar or higher metallicities should reach the Wolf-Rayet phase with small rotational velocities 
even if their 0-star progenitors rotate with speeds of up to 300 km/s (Meynet & Maeder 
2003). As massive stars age and lose their mass, most of their angular momentum is also 
carried away and stars with the strongest winds (i.e. in high metallicity environments) 
should be most affected by this spin-down process. Support for this prediction was provided 
by Harries et al. (1998) in their spectropolarimetric survey of a complete sample of northern 
WR stars. These authors found no evidence for a high frequency of highly asymmetric 
envelopes around WR stars, as would be the case for fast rotating stars (Maeder 1999). 
Their conclusions were based on the rarity of emission-line depolarization in the spectra of 
the WR stars surveyed. In fact, they found only 3 out of 16 stars showing evidence of line 
depolarization at the la level. 

In this paper, we report on our high-precision linear polarimetric observations of the 
Galactic WR binary CX Cep (WR 151) comprising an 05V star and a nitrogen-rich WN5 
star (Lewis et al. 1993). CX Cep has the second shortest period of all known WR+0 systems 
in our Galaxy (i.e. P = 2.12d) with very shallow photometric eclipses suggesting a relatively 
small orbital inclination (i.e. Lipunova & Cherepashchuk 1982 found i ~ 50" based on the 
light-curve solution). Polarimetric observations by Kartasheva (2002b) and Schulte-Ladbeck 
& van der Hucht (1989) favor, however, a much higher inchnation (i.e. i ~ 70° — 80" based 
on our reanalysis of their data; see Section 3). Kartasheva (2002b) also noted large epoch- 
to-epoch fluctuations in both the degree of matter concentration towards the orbital plane 
and the polarimetric amphtude, similar to those observed for the 1.64d WR binary CQ Cep 
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(Villar-Sbaffi et al. 2005 and Kartasheva et al. 2000). 

Our goals for this paper are to obtain a clearer picture of the matter distribution sur- 
rounding CX Cep using high-precision polarimetry with the highest time-resolution possible. 
In Section 2, we present our observations and method of reduction followed by our broad- 
band and multi-band results in sections 3 and 4, respectively. In section 5, we discuss the 
results and draw our final conclusions. 



2. Observations and Reduction 

We used the 2.1m telescope of the McDonald observatory equipped with a rotating 
Glan-prism polarimeter for our broadband multi-epoch observations (see Breger 1979). To 
maximize the photon count-rate and because wavelength-independent Thomson scattering 
is the main polarizing mechanism in the winds from massive stars, our observations were 
carried out without a filter using a Hamamatsu R943A-02 photomultiplier with extended 
red response. Each integration lasted approximately 30 minutes and included observations 
of the moonlit sky to subtract background polarization. We chose our integration times in 
order to maintain the accuracy based on photon statistics^ close to 0.03% in both Stokes 
parameters Q and U. 

Our five-band observations were carried-out at the 2.5m Nordic Optical Telescope in 
the Canary Islands using the TURPOL photopolarimeter (Piirola 1988) equipped with a 
rotating achromatic half-wave plate for linear polarimetry. This instrument allowed us to 
obtain simultaneous polarimetric observations in the UBVRI bands by using a combination 
of filters and dichroics resulting in effective wavelengths of 360nm, 440nm, 530nm, 690nm 
and 830nm for each band, respectively. The integration times were chosen such that the 
accuracy based on the least square-fit to the eight positions of the half-wave plate was less 
than 0.06% in the B band and the total integration time was less than 60 minutes. Sky 
background polarization was automatically subtracted using a plane-parallel calcite plate as 
the polarizing beam splitter (Piirola 1988). 

A summary of our observations is presented in Table 1. 

For both sites we chose the standard polarized star HD 204827 (Hsu & Breger 1982) 
to calibrate the position angle. This well-documented standard remained constant, within 



^The accuracy based on photon statistics alone (i.e. Poisson statistics) should be regarded as a lower 
limit. Bastien (1982) showed that a more accurate estimate of the actual error can be found by multiplying 
this value by 1.5 in accordance with the spurious variability of standard polarized stars. 



-4- 



~ 1°, throughout our runs and provided a rehable cahbration. For our McDonald observa- 
tions, HD 154345 (Gehrels 1974), HD 9407 (Gliese 1969) and HD 212311 (Schmidt, Elston, 
& Lupie 1992) were used as standard unpolarized stars to determine the instrumental po- 
larization. Numerous observations each night allowed us to determine an average value for 
the instrumental polarization in the Stokes parameters Q and U, which was then subtracted 
from our data. 

According to the operating manual of the TURPOL polarimeter at the NOT, the in- 
strumental polarization attributable to the optics can be considered negligible in all bands. 
This was confirmed to an accuracy (based on photon statistics) of ~ 0.03% by our own ob- 
servations of the unpolarized standards HD 154345 (Gehrels 1974), HD 9407 (Gliese 1969) 
and HD 67228 (Piirola 1977). Therefore, we did not subtract any instrumental polarization 
from our NOT observations. The details of all these cahbration measurements are presented 
in Table 2. 

Our observations were analyzed using the standard Brown, McLean, & Emslie (1978; 
hereafter BME78) Fourier method, which assumes an optically-thin, corotating matter dis- 
tribution with point sources, in order to determine the orbital parameters and the shape 
(i.e. moments) of the matter distribution. Accordingly, we represented the phase-locked 
variability of the Stokes parameters Q and U by fitting to our data (using minimization) 
a Fourier expansion up to second harmonics: 

U (0) = Uo + Ui cos A + [72 sin A + U3 cos 2A + U4 sin 2A, 
Q {4>) = Qo + Qi cos A -h (^2 sin A -h Q3 cos 2A -|- Q4 sin 2A. 

Where A = 27r0 and the light-curve phase (p was calculated using the ephemeris for 
primary minimum (i.e. WR star in front at phase 0.0) of Kurochkin (1985): 

HJDrain = 2444451.423 + 2.126897^. 

Although we expect WR binaries like CX Cep to increase their orbital period as they 
lose mass through stellar winds (see for instance the case of V444 Cyg in St-Louis et al. 
1993), P in these systems is usually less than O.lsec/yr and can be neglected. 

From the coefficients of this Fourier expansion, we calculated the orbital inclination 
i and the orientation of the line of nodes Q according to the BME78 prescriptions. Our 
definition for the orientation of the line of nodes comes from Harries & Howarth (1996) in 
which fl represents the angle of the ascending node measured from North through East with 
the constraint Q < 180° (i.e. an ambiguity of 180° exists in the determination of Q since 
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it is impossible from polarimetry to discern the ascending from the descending node). This 
definition differs from the one used by other authors, where VIqu represents the orientation 
of the major axis of the elhpse traced in the QU plane. Finally, the 7^ moments of the 
matter distribution (see BME78) allowed us to calculate TqC and t^H ^ describing respectively 
the degree of asymmetry about the orbital plane and the effective concentration of matter 
towards the orbital plane. 

The errors on the parameters were determined using a Monte- Carlo method by adding 
Gaussian noise to our data using the Poisson error associated with each point multiplied by a 
factor 1.5^ as the standard deviation of the distribution. In this manner, we generated 1000 
synthetic data sets and evaluated the orbital parameters for each set using the same Fourier 
method described above. We defined our error on a particular parameter as the interval for 
which 95% (i.e. 2a ) of our synthetic data sets could be found. 

Because CX Cep has a small orbital separation (i.e. a~ 25/?©; Lewis et al. 1993) 
and shows photometric eclipses, the system is expected to depart from the point-source 
approximation of BME78. Fox (1994) showed that a complete model for the polarimetric 
eclipses must take into account the reduction of direct unpolarized light during eclipses, the 
occultation of the shadowed regions behind the stars and the non-uniform illumination of 
electrons located close to the stars. Drissen et al. (1986) developed a first-order correction 
for the reduction of the direct unpolarized light during eclipses that they successfully used 
to correct for the phase-locked light-curve variability in the short-period WR binary CQ 
Cep (sec also Villar-Sbaffi et al. 2005). However, because CX Cep's light-curve shows only 
small fluctuations in intensity (i.e. 0.13 and 0.04 magnitude dips at primary and secondary 
eclipses respectively), we neglected this correction in our analysis. A detailed discussion on 
the departures from the BME78 point-source approximation will be presented later in this 
paper. 



3. The Polarimetric Orbit 

The goals of our broadband (Texas) observations of CX Cep were to obtain the most 
precise polarimetric orbit possible and study the epoch-to-epoch and non-BME78 polarimet- 
ric variability of the system. Unfortunately, highly variable atmospheric conditions coupled 
with the presence of the moon on our last two missions (i.e. 2001 and 2002) affected the 
reliability of our polarimetric measurements of this relatively faint system (V=12.5 mag). 
This resulted in noisier, less-reliable, data sets with fewer measurements than our 2000 data. 



^We multiplied our Poisson errors by a factor 1.5 in order to comply with the findings of Bastien (1982). 
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Our observations as a function of light-curve phase in the Stokes Q and U representation are 
shown in Figure 1. 

3.1. BME78 Analysis 

Using the BME78 diagnostics, we obtained the orbital parameters of CX Cep. The 
coefficients of our fit are presented in Table 3 along with the resulting orbital parameters. 
We also revisited the data of Kartasheva (2002b) and Schulte-Ladbeck & van der Hucht 
(1989) and applied the same method of reduction that we used for our own data. The 
errors on the polarimetric measurements by Kartasheva (2002b) are quoted as being in the 
range 0.11% to 0.15% with no information regarding the errors on each single measurement. 
Therefore, for the purpose of this analysis, we assumed an error of 0.13% on each individual 
measurement. For Schulte-Ladbeck & van der Hucht (1989), we used the reprocessed data 
(corrected using the ephemeris of Kurochkin 1985 and averaged over the two channels) 
presented in Kartasheva (2002a) and the original Poisson errors. The resulting parameters 
are also presented in Table 3. 

We note that these revised results are considerably different from those published orig- 
inally. Kartasheva (2002a and 2002b) probably used a different set of errors (i.e. weights) 
from the ones assumed in this paper to perform the minimization. This strong insta- 
bility in the solutions is reflected by the high errors on the resulting orbital parameters. 
We also note that the large differences in the Qo and Uq constants between the Kartasheva 
(2002b) V-band observations and the other unfiltered (i.e. white-light) observations is easily 
explained by the wavelength-dependence of the high interstellar polarization component (see 
Section 4) which usually peaks in the V-band. 

Our multi-epoch results for the orbital parameters agree (within 2a ) with each other 
and with the (revised) results from other authors, although our last two years are plagued by 
higher uncertainties due to the noisier data sets. Because the system had a more symmetric 
matter distribution about the orbital plane compared to other epochs (i.e. TqG <C ToH), we 
believe that our 2000 data provide a better representation of the orbit and matter distribution 
around CX Cep within the BME78 assumptions. 

3.2. Phase-Locked Non-BME78 Variability 

Phase-locked non-BME78 variability in polarimetric data has already been found in a 
few WR binaries (e.g. St-Louis et al. 1993 and Villar-Sbaffi et al. 2005). In short-period 
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binaries, the most important cause of BME78 departures is the violation of the point-source 
assumption due to the comparable size of the orbital separation and the stellar radii (e.g. 
a ~ 25i?0 and Rq + RwR ~ IS-R© for CX Cep). 

In order to better appreciate these deviations from the predicted double-wave behavior 
in CX Cep, we binned our 2000 data (i.e. the most reliable data set) in phase intervals of 
0.03. This binning allowed us to average-out the stochastic polarimetric fluctuations caused 
by the presence of fast-moving blobs in the wind of the WR star (see Robert et al. 1989)'^. 
The effect of these blobs on the polarization should be maximal around phases 0.0 and 0.5 
due to the alignment along the line-of-sight of the stars during eclipses. This is in agreement 
with the larger scatter observed in our data around those phases. Finally, we rotated the 
data by —^qu — 30° in the QU plane^ to the orbital plane of symmetry of the system and 
plotted the polarization in Stokes Q and U as a function of light-curve phase (Figure 2). 
A clear non-BME78 (i.e. non-double-wave) structure appears around phase 0.5 when the 
O star eclipses the WR star. A similar structure has been observed in the WR binaries 
V444 Cyg (St-Louis et al. 1993) and CQ Cep (Villar-Sbafh et al. 2005). This non-BME78 
modulation was successfully modeled for V444 Cyg by considering the effect of the gradual 
eclipse by the O star of the dense inner regions of the WR wind during secondary eclipse. 
This simple model allowed the authors to constrain the stellar radii of both stars along with 
the mass-loss rate and velocity law of the WR star. 

Following this success, we decided to develop a similar model to simulate the polarimetric 
variability in CX Cep not only at the eclipses but over the entire double- wave orbit. Our 
main goal is to determine the mass-loss rate of the WR star. We used the same formalism 
as St-Louis et al. (1993) and solved the modified BME78 integrals that assume an optically- 
thin wind corotating in the frame of the binary. The finite size of the stars was accounted for 
by setting the electron density equal to zero in the shadowed regions behind the stars and 
correcting for the depolarization of photons scattered close to the stars using the Cassinelli 
et al. (1987) factor. 

The orbital separation of 25i?0 from Lewis et al. (1993) was used and the magnitude 
difference M^r — Mq = 1.4 was taken from van der Hucht (2001). We assumed that the 
WR wind consists exclusively of ionized helium and followed the prescriptions of Moffat & 
Marchenko (1993) to calculate the number of free electrons per nucleon in the wind for a 



■^Here we assume that the spatial distribution of blobs over a long period of time (i.e. 2 weeks) follows 

the density distribution of the wind. 

''This value was chosen from our most reliable 2000 data and is related to the orientation for the line of 
nodes O by the relation VIqu = 20 — 180. 
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typical WNE star. The radius of the O star was kept fixed at 11 Rq based on the most 
recent models of Martins et al. (2005) for an 05V star. Although the spectral type of the 
O star companion is still a source of debate (Kartasheva 2002a), our models were not very 
sensitive to changes in the 0-star radius of ±2Rq (see Table 4). 

We initially fixed the orbital inchnation of our models at 73" based on the 2000 broad- 
band results. However, wc quickly realized that a lower inclination produced a much better 
fit to the binned data. An inclination of 65° was finally adopted. This value differs by Aa 
from the one determined by the BME78 analysis. However, this discrepancy agrees with the 
results of Fox (1994) who showed that by considering occultation effects in the envelope of 
a modeled binary system, deviations of ~ 3" could be found in the inclination derived by 
Fourier analysis of a synthetic polarimetric orbit. 

The integrals were solved within a volume equal to 10 times the orbital separation by 
Monte-Carlo integration using a Sobol quasi-random number sequence in three dimensions 
with a word-length of 30 bits. The major improvement of our model was to replace the 
cut-off radius used by St-Louis et al. (1993) by the optical-depth r along the line-of sight, 
therefore replacing an arbitrary parameter (i.e. the cut-off radius) by a physical one. The 
cut-off radius was initially introduced in order to exclude high density regions in the inner 
wind where multiple-scattering leads to significant depolarization. 

The optical-depth for a scattered photon at coordinates {xo, Vo, Zo) was calculated by 
assuming that the only source of opacity was Thomson scattering: 



Where ax — 6.65 x 10~^^cm~^ is the Thomson scattering cross-section. The electron 
density ne{x,y,z) was determined using the principle of mass conservation at every point 
in the accelerated wind using a /3-type velocity law with (3 — 5 and a sonic radius (i.e. 

Rs) of 2Rq for a typical WNE star. The high value of (3 and the small sonic radius of 
the wind follow the optically-thick radiation-driven wind models of Nugis & Lamers (2002) 
who were able to reproduce the high mass-loss rates of Wolf-Rayet stars^. Tmax = 2/3 was 
chosen as the maximum optical-depth beyond which a photon became unpolarizcd and did 
not contribute to the total polarization. This optical depth corresponds to the adopted 
definition of a star's photosphere and resulted in a good fit to the data (see models MIO and 



lower /J-value with a larger sonic radius (i.e. (3 = 1 and Rg = 4:Rq) was also found to fit our data 
without affecting considerably the determined mass-loss rate. 



oo 
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Mll in Table 4). However, because we assume that Thomson scattering is the only source of 
opacity and neglect line transitions, we are in fact underestimating the photospheric radius 
and overestimating its Thomson scattering optical-depth (i.e. Tmax)- As our models will 
prove (see Table 4), a lower T^ax will result in a higher mass-loss rate for the WR star. 
Therefore, the mass-loss rate determined with Tmax = 2/3 represent a lower limit compared 
to the actual mass-loss rate. 

In Figure 2 (solid line) and Table 4 (model Ml) we present the best fitting model to our 
binned data. The model was fitted by minimization by first determining the horizontal 
(i.e. A0) and vertical (i.e. AQ and AU) off-sets and finally by fitting the value 7 = Af /\4o 
which controls simultaneously the ampUtudes in Stokes Q and U (see St-Louis ct al. 1993). 
Although our best fit is quite good, there are still a few structures in Stokes Q at phases 
0.0 and 0.5 that could not be reproduced. These structures are probably signs of a complex 
matter distribution which might include a distorted wind-wind collision zone and a disk-like 
density enhancement of the wind (see section 4.2). To avoid any problems caused by these 
structures on our final results we excluded from the minimization any data points in 
Stokes Q located within ±0.1 in phase around the eclipses. 

Assuming a characteristic terminal velocity for a WN5 wind of 1700 km/s, we obtain 
a mass-loss rate for the WR component of 0.3 x 10~^ Mq/ut. A few other models showing 
the effect of varying the orbital inclination (models M2 and M3), the exponent f3 of the 
velocity law (models M4 and M5), the star's radii (models M6, M7, M8 and M9) and the 
maximum optical depth (models MIO and Mil) are also shown in Table 4. Surprisingly, the 
mass-loss rate determined from all the models is remarkably stable and only fluctuates by 
about 0.2 X 10~^ Mq/uv around the value determined from our best fitting model. As we 
stated above, the adopted maximum Thomson scattering optical depth (i.e. Tmax — 2/3) is 
our main source of uncertainty and, in fact, results in a lower limit for the mass-loss rate 
of the WR star. However, because spectral lines only amount to approximately 10% of the 
total flux from a WR star, we do not believe that Tmax could be considerably lower than our 
assumed value of 2/3. 

We also note that a small bias towards low mass-loss rates is inherent to our polarimetric 
model since we neglect the weaker 0-star wind. In relation to WR stars, 0-stars have lower 
mass-loss rates and faster winds and, therefore, lower densities in their circumstellar matter 
distribution. Since the fltted parameter 7 (i.e. 7 = M/Voo) from which we derive the mass- 
loss rate is actually an average over the whole wind, the low-7 contribution from the 0-star 
wind will decrease the mass-loss rate determined from our model. Fortunately, the matter 
distribution surrounding CX Cep is dominated by the strong WR wind (see Lewis et al. 
1993, their Figure 14) and this bias is probably not very high. 
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Finally, based on the resulting mass-loss rates from all our models, we believe that the 
mass-loss rate of the WR star is in the range 0.3 — 0.5 x lO~^M0/|/r, with a preference 
towards the higher value. A detailed two-wind model including wind-wind collision effects 
would certainly provide a more accurate result. 



By measuring the constant intrinsic polarization (CIP) of a star, we are in fact measuring 
the degree of spherical asymmetry of the matter distribution surrounding it. A high CIP 
implies a strong departure from spherical symmetry. In order to determine the CIP of a star 
we must fit and remove the contribution of the interstellar polarization to the total observed 
polarization. For massive stars, where (wavelength-independent) Thomson scattering is the 
dominant source of linear polarization, this is usually done by obtaining multi-wavelength 
(i.e. UBVRI) polarimetric observations of the star in order to subtract the wavelength- 
dependent interstellar contribution. Because in this section we are only interested in the CIP 
of the system, our multi-band observations are scarcer in phase coverage and serve essentially 
to determine the zero-point of the polarimetric variability in each band (see Figure 3) . 

We note that a small wavelength-dependent intrinsic contribution to the polarization 
is expected due to the presence of strong emission fines in the spectrum of a WR star (see 
for instance Moffat & Piirola 1993). These lines form far from the dense central regions 
of the wind and their polarization can differ from that of the surrounding continuum (see 
Villar-Sbaffi et al. 2005). However, tfie line-fiux in a standard Jofinson band only represents 
~ 10% of tfie total flux from tfie star and we decided to neglect tfiis effect in our analysis. 
Tfie results of tfie BME78 analysis of tfie multi-band observations are presented in Table 5. 

Tfie CIP of CX Cep was determined using tfie metfiod described in Villar-Sbaffi ct al. 
(2005). We assumed tfiat tfie constant interstellar polarization could be represented by a 
modified Serkowski's law (Serkowski, Matfiewson, & Ford 1975): 



Wfiere V^ax is tfie maximum interstellar polarization found at wavelengtfi Xmax^ K (Amax) = 
0.01 + 1.66Ama^ (Wfiittet et al. 1992) and Ois (A) = Ci + C2/A (Dolan & Tapia 1986, Moffat 
& Piirola 1993 and Matsumura et al. 2003). Tfie constant Stokes parameters Qo and Uq 
from tfie BME78 analysis can tfien be expressed as: 



4. 



The Constant Intrinsic Polarization 



Uis{\) = PmaxSinI 

Qls{X) = PmaxCOS 



[2ejs] e-^(^— 
:[2eis] e-^(^— ^^'^'l^). 
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1+cos-^ i 
1+cos-^ i 



- + Ujs (A) , 
i + Qis (A) . 



Where Uwr and Qwr are the components of the CIP vector and C/3, Q3 and i were 
taken from our fit to the broadband 2000 data (see Table 3). 

Using these equations and a non-linear Levenberg-Marquardt code, we fitted simulta- 
neously our values of Qo and Ug to obtain the CIP of the system (i.e. Qwr and Uwr) and 
the parameters describing the interstellar polarization (i.e. Pmax, ^max, Ci and C2). The 
numerical weight of each band was taken to be the rms error on the double- wave fit to Q(0) 
and U(0) found in Table 5 and not the actual errors on Qo and Uo determined from the fit. 
We justify this statement by the fact that the constants Qo and Uq only represent the sum of 
the interstellar and intrinsic polarizations to the extent that the BME78 assumptions hold 
for this system. Therefore, we believe that the rms errors of the double-wave fit represent 
more accurately the relative weight to be given to each band. The resulting parameters for 
this first model (i.e. Model 1) are presented in Table 6 and the fitted curve along with our 
data can be found in Figure 4. 

Because of the possible instability of our solution due to the weak A-dependency of the 
interstellar polarization, we performed another fit with the constraint Qwr = Uwr = 0% to 
test the possibility that our first solution was statistically insignificant. The parameters of 
this second fit and the corresponding curve (i.e. Model 2) can also be found in Table 6 and 
Figure 4. 

With an rms error less than half that of Model 2, Model 1 provides a much better fit 
to the constant polarization of CX Cep^. We can therefore say with certainty that the CIP 
in not zero and probably very high. However, the resulting interstellar polarization seems 
very high. According to Serkowski, Mathewson & Ford (1975), the maximum interstellar 
polarization observed for completely aligned dust grains in a purely regular external magnetic 
field is related to the extinction E(B-V) along the line-of-sight by the equation Pmax — 
9E{B — V). Using the absorption Ay = 3.77 mag in van der Hucht (2001) and the commonly 
accepted ratio Ry = 3.1 between the extinction and absorption in the V band we obtain 



^The rms error presented in this paper is calculated using the formula: rms 



VYl {Vi ~ yfit{xi)Y /{N — M), where M is the number of fitted parameters. This rms error gives an 
i=l 

unbiased measure for the quality of the fit, which takes into account the reduced number of free parameters 
in Model 1 compared to Model 2. 
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a maximum interstellar polarization along the line-of-sight towards CX Cep of almost 11%. 
This is larger than the value determined from Model 1 and, therefore, in accordance with 
our results although such a high polarization would necessitate an extreme case of dust-grain 
alignment. Another source of supporting evidence for this high interstellar polarization is 
the presence of a normal (i.e. should not show intrinsic polarization) B0.5V star within 
0.11" of CX Cep at a distance of 2.3kpc (cf. ~ 6kpc for CX Cep; van der Hucht 1988) with 
an average interstellar polarization of 4.5% oriented at an angle of 41" (Heiles 2000). This 
angle is the same as that determined from Model 1 but differs by 10° from the interstellar 
polarization angle of Model 2, therefore supporting the validity of Model 1. 

The main source of uncertainty on our CIP determination is probably the scarce data 
coverage centered mainly on phases 0.0 and 0.5 where the scatter is larger (see Figure 1) due 
to the presence of blobs in the wind and the alignment of the light-sources during eclipses. 
However, since our multi-band observations cover more than three orbits (i.e. 8 nights) and 
the eclipses were therefore observed at least three times, we believe that any epoch-to-epoch 
fluctuations were averaged-out (i.e. assuming again that the blob distribution follows that 
of the wind over a long period of time). An estimate of the error caused by the stochastic 
polarimetric fluctuations can be obtained through the covariance matrix produced by the 
Levenberg-Marquardt method used to fit and remove the interstellar polarization. The 
diagonal elements of this matrix are the squared uncertainties on the fitted parameters. 
These uncertainties depend on the weights (i.e. errors) given to the data points that we 
assumed to be the rms errors on the BME78 double-wave fits. These rms errors measure 
the scatter in Q and U caused by blobs (and other non-BME78 behavior) around the mean 
BME78 curve; thus, they provide an estimate of the maximum vertical shift (i.e. in Stokes 
Q and U) that blobs can impart to our data. Therefore, the covariance matrix provides an 
upper-limit estimate of the error on our fitted CIP. Using this method we obtain a maximum 
error on the intrinsic polarization of Upj^^ = 1%, which is almost four times lower than the 
observed CIP (i.e. P = 3.89%). 

Based on these facts, it appears that CX Cep has a very high intrinsic polarization (i.e. 
PwR ~ 3 — 4%) aligned at 26° from the equatorial plane of reference. This angle doesn't 
correspond to either the orientation of the orbital plane (i.e. Q) or the polar axis (i.e. O-I-90) 
of the system. A similar phenomenon was recently observed in the short-period WR binary 
CQ Cep (Villar-Sbaffi et al. 2005) and indicates an apparent lack of correlation between 
binarity and wind geometry in WR-l-0 systems. To confirm this assertion, we performed a 
fit on the constants Qo and Uo with the constraint 9wr — 0, + 90", as would be the case for 
a tidally distorted wind. The results of this third model have an rms error more than twice 
that of model 1 and confirm the misalignment of the CIP in relation to the orbital axis (see 
Table 6). 
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5. 



Discussion and Conclusions 



Although it has the second shortest period of all know WR+0 binaries in our Galaxy, 
CX Cep is a relatively poorly studied system. In this paper, we have confirmed the po- 
larimetric orbital parameters determined by other authors and considerably improved their 
accuracy. The BME78-derived orbital inchnation confirmed the almost edge-on configuration 
of the system and provided an accurate orientation for the line of nodes [i.e. i — (73 ± 2)° 
and Q — {75 ± 1)°]. However, by including the finite stellar size effects in a model of the 
polarimetric orbit, we found a lower and more accurate inclination of 65°. Assuming an error 
on the orbital inclination of 3°, the resulting masses for both components can then be found 
using the orbital solution of Lewis et al. (1993): 



The mass of the 0-star agrees with the spectroscopic masses for Galactic 05-6V stars 
(Martins et al. 2005) and the WR mass falls within the acceptable range for a rotating WNE 
star with an initial mass of 60i?o (Meynet & Maeder 2003). 

Our polarimetric model also allowed us to obtain the mass-loss rate of the WR com- 
ponent, which was found to lie in the range 0.3 — 0.5 x W^MQ/yr and was remarkably 
independent of the adopted input parameters. However, this value is quite low compared 
to other WN5 stars with mass-loss rates close to 1 x 10~^MQ/yr. An earlier (i.e. hotter) 
spectral type (i.e. WN3-4) would be more appropriate for such a low M (Nugis & Lamers 
2002). In fact, Lewis et al. (1993) concluded that a WN4 spectral type was also acceptable 
for the WR component in GX Gep. 

Our most surprising discovery is the presence of a large 3-4% GIP in GX Gep, higher 
than that of the most polarized Be stars (Poeckert et al. 1979). Although the uncertainty 
on this value (and the corresponding high interstellar polarization) is probably high (i.e. up 
to 1% in the worst case scenario), our analysis shows that the probability that the wind is 
spherically symmetric is very low. Based on the radiative transfer model presented in Harries 
et al. (1998), this GIP imphes an equator-to-pole density ratio higher than 5. Like for the 
WR binary GQ Gep (Villar-Sbaffi et al. 2005), this departure from spherical symmetry is 
characterized by a GIP vector misaligned in relation to the natural axes of symmetry of a 
binary (i.e. the polar and equatorial axes) which points towards a rotational origin for the 
asymmetry. Alternatively, this could reveal the presence of a magnetic field which is affecting 
the symmetry of the wind. 

Within the observations of Harries et al. (1998) and the models of Meynet & Maeder 



Mo = 

MwR 




sin i 



© 
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(2003), CX Cep is an anomaly. Massive stars at solar metallicity should reach the WR phase 
with low rotational velocities and, therefore, spherically symmetric winds. For this reason, 
the extremely flattened wind of CX Cep should provide valuable information regarding the 
effects of rotation on massive star evolution and should be studied further. 

AVS, NSL and AFJM are grateful for financial support from NSERC (Canada) and 
FQRNT (Quebec). We would also like to thank Andrei Berdyugin for his help operating the 
NOT telescope and for sharing his observing time with us. 
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Fig. 1. — Normalized Q and U Stokes parameters for our unfiltered McDonald observations 
as a function of the 2.12d light-curve phase of CX Cep. These observations were fitted 
separately for each epoch to a Fourier expansion up to second harmonics (solid line) following 
the prescriptions of BME78. 
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Fig. 2. — Binned broadband observations of CX Cep as a function of light-curve phase 
rotated by —^Iqu to the orbital plane of symmetry. The solid line represents our best fitting 
model (i.e. model Ml from Table 4) accounting for the polarimetric eclipse of the WR wind. 
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Fig. 3. — Normalized Q and U Stokes parameters from our multi-band (UBVRI) NOT 
observations versus the 2.1 2d light-curve phase of CX Cep. These observations were fitted 
to a Fourier expansion up to second harmonics (solid line) following the prescriptions of 
BME78. 
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Fig. 4. — Constant components of the polarization Qo and Uo as a function wavelength 
obtained by fitting our NOT observations with a Fourier function up to second harmonics. 
The solid line represents our fit using a modified Serkowski function plus an additional 
wavelength-independent term corresponding to the intrinsic polarization of the system. The 
dashed line represents the same fit obtained with the constraint Qwr — Uwr — 0% (i.e. for 
a spherically symmetric wind). 
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Table 1. Summary of Observations 





Epoch 1 (2000) 


Epoch 2 (2001) 


Epoch 3 (2002) 


Observatory 


McDonald 


NOT 


McDonald 


McDonald 


Telescope size 


2.1m 


2.5m 


2.1m 


2.1m 


Instrument 


Breger polarimeter 


TURPOL 


Breger polarimeter 


Breger polarimeter 


Filters 


None 


UBVRI 


None 


None 


First night (UT) 


August 31"* 


October 8*'' 


August 30*'' 


September lO"' 


Number of nights 


14 


20 


14 


14 
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Table 2. Calibration Measurements 



McDonald 




NOT (2000) 




2000 2001 2002 


U 


B V R 


I 


P^nstrumental'' 0.081% 0.083% 0.084% 

- OpubHshed'' 42.7° 41.8° 40.4° 


0.029% 
-71.3° 


0.018% 0.023% 0.021% 
-71.7° -71.7° -71.9° 


0.069% 
-70.6° 



^These values represent the average observed polarization of the assumed unpolarized stars HD 

154345 (Gchrels 1974), HD 9407 (Gliese 1969) and HD 212311 (Schmidt, Elston, & Liipic 1992) 
for our McDonald (unfiltered) observations and HD 154345 (Gehrels 1974), HD 9407 (Gliese 1969) 
and HD 67228 (Piirola 1977) for our NOT observations. 

'^HD 204827 was chosen as our standard polarized star for both sites with OpuUished = 59.3° 
according to Hsu & Breger 1982. 



Table 3. Polarimetric Parameters of CX Cep from our Unfiltered Observations 







Unfiltered Results (McDonald) 




Schultc-Ladbcck & 
van der Hucht (1989) 
Unfiltered 


Kartashcva 
(2002b) 
V band 


2000 


2001 


2002 


Combined 


Uo \/0) 


D.UDo it U.UUz 


6.176 ± 0.004 


6.076 ± 0.007 


(\ not; 4- n nni 
D.uyo It U.UUi 


(\ ^ A -\- f\ n/1 
D.14 It U.U4 


(.U4 it U.U/ 


TTi 

Ul \/o) 


n 099 -1- n nn9 

-U.UZZ III u.uuz 


-0.047 ± 0.007 


0.075 ± 0.006 


n ni -t- n nno 
-u.uio m u.uuz 


n ns -t- n nfi 

U.UO It U.UD 


n 1 1 -1- n n7 
-u.ii m u.U( 


U2 


n 09*7 -U n nn9 
-U.Uz( it U.UUz 


-0.003 ± 0.003 


0.049 ± 0.010 


-u.uuy It u.uuz 


n nfi 4- n nK 
-U.UO It U.UD 


n 1 a -J- n cxi 

-U.lcS It U.U* 


U3 [/<3) 


U.UOO it U.UUz 


0.057 ± 0.006 


-0.023 ± 0.006 


U.UlcS It U.UUz 


n nc 4- n nK 
U.UcS It U.UD 


U.UD it U.U4 


U4 (%) 


0.182 ± 0.003 


0.196 ± 0.005 


0.185 ± 0.007 


0.188 ± 0.002 


0.17 it 0.06 


0.31 it 0.09 


Qo (%) 


-0.390 ± 0.002 


-0.400 ± 0.004 


-0.476 ± 0.007 


-0.413 it 0.001 


-0.01 it 0.04 


-0.08 it 0.07 


Qi (%) 


0.006 ± 0.002 


-0.162 ± 0.007 


-0.034 ± 0.006 


-0.025 it 0.002 


-0.05 it 0.05 


0.13 it 0.07 


Q2 (%) 


-0.006 ± 0.002 


0.047 ± 0.003 


-0.046 ± 0.010 


-0.003 it 0.002 


-0.04 it 0.06 


0.04 it 0.07 


Q3 (%) 


-0.240 ± 0.002 


-0.180 ± 0.006 


-0.153 ± 0.006 


-0.200 it 0.002 


-0.12 it 0.05 


-0.25 it 0.04 


Q4 (%) 


-0.056 ± 0.003 


-0.119 ± 0.005 


-0.062 ± 0.007 


-0.076 it 0.002 


0.04 it 0.05 


-0.19 it 0.09 




73±2 


79±3 


63±20 


70it2 


68itl8 


79itl2 


Q}" (°) 


75±1 


69±2 


61±6 


72itl 


50it38 


67itl3 


ToG'' (%) 


0.018 


0.087 


0.053 


0.015 


0.067 


0.127 


ToH'' (%) 


0.250 


0.269 


0.165 


0.219 


0.166 


0.397 


rms[/ (%) 


0.059 


0.071 


0.065 


0.082 


0.170 


0.150 


rmsQ (%) 


0.075 


0.080 


0.077 


0.103 


0.121 


0.177 



''The errors on TqG and ToH were found to be less than 0.02% at the 2a level and we decided not to present them explicitly. 
''Our definition for f2 is that of Harries & Howarth (1996). 



Table 4. Parameters of the Non-BME78 models of CX Cep 





Ml 


M2 


M3 


M4 


M5 


M6 


M7 


M8 


M9 


MIO 


Mil 


in 


65 


60 


70 


65 


65 


65 


65 


65 


65 


65 


65 


13 


5 


5 


5 


4 


6 


5 


5 


5 


5 


5 


5 


Ro (i?©) 


11 


11 


11 


11 


11 


9 


13 


11 


11 


11 


11 


Rs (Rq) 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


1.0 


3.0 


2.0 


2.0 




0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.33 


1.00 


A4) 


0.033 


0.033 


0.033 


0.033 


0.032 


0.032 


0.035 


0.033 


0.031 


0.031 


0.033 


AU (%) 


4.96 


4.96 


4.95 


4.96 


4.95 


4.96 


4.95 


4.97 


4.95 


4.95 


4.95 


AQ (%) 


-3.74 


-3.70 


-3.77 


-3.74 


-3.75 


-3.73 


-3.74 


-3.76 


-3.72 


-3.72 


-3.73 


7(10-"MQ/A;m) 


6.1 


5.4 


6.6 


6.6 


5.6 


5.5 


6.9 


8.4 


4.9 


9.3 


4.8 


M(10-5AfQ/j/r)^ 


0.33 


0.29 


0.36 


0.36 


0.30 


0.30 


0.37 


0.45 


0.27 


0.50 


0.26 


rms 


0.040 


0.043 


0.047 


0.042 


0.043 


0.043 


0.043 


0.051 


0.053 


0.052 


0.043 



^The mass-loss rate was determined assuming a typical terminal velocity (i.e. Voo) for the wind of a WN5 star 

of 1700 km/s. 

''This rms represents the average scatter of our binned data points from the modeled curve calculated by 
excluding any points in Stokes Q located at ±0.1 in phase from the eclipses . 



Table 5. Polarimetric Parameters of CX Cep from our Multi-Band Observations 





U 


B 


V 


R 


I 


Uo (,/c) 


0.4o ± U.U4 


D.4o it U.U4 


D.oi it U.U4 


a AT -L n nQ 
D.4 ( it U.Uo 


/1 1 J- n nQ 
0.41 it U.Uo 


Ul (/o) 


n ri/i -I- n no 


n nv -U n no 
U.U 1 it U.UZ 


n nv -U n no 
U.Ui it U.U^ 


n HQ _i_ n ni 
-U.UvS it U.UI 


n no _i_ n ni 
-U.yjZ it U.UI 


TT„ (Of") 




u.uy It U.U4I: 


U.UD It U.UO 


n 1 1 -1- n n'^ 

U.±l It U.Uo 


n 1 Q _|_ n n'^ 
-U.lo It U.UO 


U3 (/o) 


u.iy it u.uo 


U.Uo it U.U4 


n n? -U n nc; 
-U.U/ it U.UO 


nil -U n m 
U.ll it U.Uo 


n OQ -U n r\A 

U.^o it U.U4 


U4 (%) 


0.17 ± 0.03 


0.19 ± 0.03 


0.23 ± 0.04 


0.09 ± 0.02 


0.02 ± 0.02 


Qo (%) 


-0.29 ± 0.04 


-0.18 ± 0.03 


-0.42 ± 0.03 


-0.47 ± 0.03 


-0.75 ± 0.03 


Qi (%) 


0.03 ± 0.02 


-0.02 ± 0.02 


-0.01 ± 0.01 


-0.09 ± 0.01 


-0.05 ± 0.01 


Q2 (%) 


0.16 ± 0.04 


-0.07 ± 0.04 


0.09 ± 0.04 


0.01 ± 0.03 


0.03 ± 0.03 


Q3 (%) 


-0.18 ± 0.04 


-0.38 ± 0.03 


0.00 ± 0.04 


-0.05 ± 0.03 


-0.11 ± 0.03 


Q4 (%) 


-0.18 d= 0.03 


-0.02 ± 0.03 


-0.20 ± 0.03 


-0.21 ± 0.02 


-0.04 ± 0.02 




89±7 


76±14 


86±9 


81±8 


87±8 




67±7 


82±8 


65±7 


74±7 


58±8 


ToG^ (%) 


0.10 


0.07 


0.07 


0.08 


0.09 


ToH- (%) 


0.36 


0.37 


0.31 


0.24 


0.26 


rms(7 (%) 


0.16 


0.18 


0.17 


0.09 


0.15 


rmsQ (%) 


0.18 


0.19 


0.18 


0.06 


0.24 



to 



I 

^Tlic errors on TqG and ToH were found to be less than 0.02% at the 2<t level and we 
decided not to present them explicitly. 

''Our definition for H is that of Harries & Howarth (1996). 
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Table 6. Interstellar and Intrinsic Components of the Polarization from CX Cep 





Model 1 


Model 2 


Model 3 


'^max (^) 


10.17 


6.75 


6.09 


Xmax (nm) 


545 


557 


554 


ci n 


41 


51 


46 


C2 ("nm) 


-31 


-52 


-58 


VwR (%) 


-3.05 


0* 


0.60 


QWR (%) 


-2.36 


0* 


-1.03 


PWR (%) 


3.89 


0* 


1.19 




26 


N/A 


165* 


rms (%) 


0.09 


0.22 


0.24 



* Fixed. 



